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Abstract In this paper we review the search for astrophysical neutrinos. We begin by summa- 
rizing the various theoretical predictions which correlate the expected neutrino flux with data 
from other messengers, specifically gammas and ultra-high energy cosmic rays. We then review 
the status and results of neutrino telescopes in operation and decommissioned, the methods 
used for data analysis and background discrimination. Particular attention is devoted to the 
challenge enforced by the highly uncertain atmospheric muon and neutrino backgrounds in re- 
lation to searches of diffuse neutrino fluxes. Next, we examine the impact of existing limits on 
neutrino fluxes on studies of the chemical composition of cosmic rays. After that, we show that 
not only do neutrinos have the potential to discover astrophysical sources, but the huge statistics 
of atmospheric muons can be a powerful tool as well. We end by discussing the prospects for 
indirect detection of dark matter with neutrino telescopes. 
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1 INTRODUCTION 

It has long been recognized that high energy protons produced in cosmic ray 
accelerators would also generate an observable flux of cosmic neutrinos, mainly 
through charged pion production in collisions with the ambient gas or with radi- 
ation fields (1). Viceversa, the emission of neutrinos necessarily implies the ex- 
istence of relativistic baryons, and consequently the acceleration of cosmic rays. 
Upcoming catalogues of cosmic ray and neutrino sources would thus in principle 
be identical were it not for the different properties of the messengers themselves. 

Neutrinos can escape from the innermost regions of galaxies where light and 
other kinds of electromagnetic radiation are blocked by matter and hence can be 
tracers of processes which stay hidden to traditional photon astronomy. Further- 
more, neutrinos carry with them information about the site and circumstances 
of their production. Undeflected and unabsorbed they can point back to very 
distant sources, resolving the origin of the highest energy cosmic rays and the 
underlying acceleration mechanism. Besides, unlike gamma-rays that can be 
produced via synchrotron radiation, or inverse Compton scattering, neutrinos 
provide incontrovertible evidence for acceleration of baryonic cosmic rays. 

MeV neutrino astronomy has been possible for about 45 years (2). Thus far 
two sources of non-atmospheric neutrinos have been identified: the Sun (3) and 
SN1987a (4). During the next decade we will be able to observe the universe 
using high energy "neutrino light." This will doubtless illuminate a wondrous 
new view of the universe. 

Candidate sources of high-energy neutrinos are of both galactic and extra- 
galactic varieties, the latter being expected to dominate at the highest energies 
just as for the parent cosmic rays. Energy balance considerations suggest that 
the diffuse galactic flux of neutrinos predominantly originates in supernova rem- 
nants (5). Dense molecular clouds, often found in star forming regions where 
the supernovae (SNe) explode, are particularly efficient at converting protons 
into pions that decay into gamma rays and neutrinos (6). Other possibilities for 
galactic sources include microquasars (7), 7-ray binaries (8), and pulsar wind neb- 
ulae (9). By far, the most likely extra-galactic sources are active galactic nuclei 
(AGN) (12) and gamma-ray bursts (GRBs) (13). [For reviews see e.g., (10) and 
for neutrino flux estimates see (11)]. In addition, the so-called cosmogenic neu- 
trinos are guaranteed to exist, their amount being still highly uncertain. These 
originate from the decay of charged pions produced in the interactions of protons 
with the isotropic photon background (14-16). 

As for conventional astronomy, neutrino astronomers observe the neutrino sky 
through the atmosphere. This is a curse and a blessing; the background of neu- 
trinos produced by cosmic rays in interactions with air nuclei provides a beam 
essential for calibrating the detectors and demonstrating the neutrino measure- 
ment technique. It also provides an opportunity to probe neutrino standard 
oscillations and those arising from new physics, such as violation of Lorentz in- 
variance (17). Especially unique is the energy range of the atmospheric neutrino 
beam between 10 GeV < E u < 10 8 GeV, not within reach of accelerators, that 
may reveal the existence of neutrinos from prompt decays of charmed mesons and 
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baryons (18). Over the next decade, a data set of the order of one million atmo- 
spheric neutrinos will be collected. The statistics will be so large that mapping 
of the Earth's interior will be possible via neutrino tomography (19). 

Though neutrino astrophysics is the central motivation to build neutrino tele- 
scopes, ultrahigh energy cosmic neutrinos are also unique probes of new physics 
as their interactions are uncluttered by the strong and electromagnetic forces and 
can reach center-of-mass energies of y/s ~ 245 TeV (20). Rates for new physics 
processes, however, are difficult to test since the flux of cosmic neutrinos is virtu- 
ally unknown. It is possible to mitigate this by using multiple observables which 
allow one to decouple effects of the flux and cross section (21). Large neutrino 
telescopes may also open a new perspective for indirect searches of dark matter 
by detecting neutrinos produced via annihilation of weakly interactive massive 
particles trapped within the Sun, or at the center of the Earth (22). 

The layout of the paper is as follows: in Sec. 2 we overview the connection 
between cosmic ray observations and expected neutrino fluxes; in Sec. 3 we re- 
view the status of neutrino telescope construction, operation and data analysis; 
in Sec. 4 we analyze the upper limits on neutrino fluxes set by the different ex- 
periments to assess their scientific impact; in Sec. 5 we provide a summary and 
discuss future prospects. 

2 COSMIC RAY ^ NEUTRINO CONNECTION 

The cosmic ray flux falls as an approximate power-law in energy, Jcr(E) oc E~ a , 
with a ~ 2.7 from about 10 GeV up to the 'knee' in the spectrum at E ~ 

3 x 10 6 GeV where it steepens to a ~ 3.1; it then steepens further to a ~ 3.25 
at E ~ 4 x 10 8 GeV (the so called 'dip') and flattens back to a ~ 2.69 at 
E ~ 3 x 10 9 GeV, the so-called 'ankle' (23). The spectrum extends up to at least 
E ~ 10 11 GeV with a steep spectrum {a ~ 4.2) beyond E = 5x 10 10 GeV (24), 
in accord with the prediction by Greisen, Zatsepin and Kuz'min (GZK) (25). 

The energy density associated with the flux of galactic cosmic rays is pe ~ 
10~ 12 erg cm -3 (26). This is about the value of the corresponding energy density 
B 2 /8ir of the microgauss magnetic field in the Galaxy. Given that the average 
containment time of the cosmic rays in our Galaxy is ~ 3 x 10 6 yr, the power 
needed to maintain this energy density is ~ 10 -26 erg/cm 3 s. For a nominal 
volume of the galactic disk of 10 67 cm 3 this requires an accelerator delivering 
power of about 10 41 erg/s. This happens to be 10% of the power produced by 
supernovae, releasing 10 51 erg every 30 yr. This coincidence is the basis for the 
idea that shocks produced by supernovae expanding into the interstellar medium 
are the origin of the galactic cosmic rays. 

The conversion of 10 50 erg of energy into particle acceleration is believed to 
occur by diffusive shock acceleration in the young (10 3 — 10 4 yr) SN remnant 
expanding into the interstellar medium. If high energy cosmic rays are indeed 
associated with the remnant, they will interact with hydrogen atoms in the in- 
terstellar medium to produce pions that decay into roughly equal numbers of 
photons and neutrinos: the former via 7r° — > 77, the latter via 7r + — > e + u e u^u lx 
(and the conjugate process). These may provide us with indirect evidence for 
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cosmic ray acceleration. The observation of these pionic gamma rays has been 
one of the motivations for neutrino astronomy as well as, to a large extent, for 
gamma-ray astronomy 

Whereas the details are complex and predictions can be treacherous, a simple 
estimate of the gamma ray flux associated with a supernova remnant can be made 
following Aharonian, Drury and Volk (27). The emissivity in pionic gamma rays 
is simply proportional to the density of cosmic rays for energy in excess of 1 TeV, 
^Cr(> 1 TeV) ~ 4 x lCT 14 cm~ 3 , and to the target density of hydrogen atoms, 
nn- The proportionality factor is determined by particle physics. 

The emissivity (number /volume/time/energy) of neutral pions resulting from 
an isotropic distribution of highly relativistic particles following a power-law en- 
ergy spectrum dncR(Ejsr)/dEN oc E^ a is given by 

Q K (E n ) = cn R J d ^{E N )^{E 7 ,,E N )dE N (1) 

where E t ^(E 7r ) is the minimum energy per nucleon required to produce a pion 
with energy E n , and daA.(E n ,E^)/dE n is the differential cross section for the 
production of a pion with energy E n in the lab frame due to the collision of a 
nucleus A of energy per nucleon En with a hydrogen atom at rest. The differential 
cross-section can be parametrized by 

^(E n , E N ) ~ ^ xF n (x, E N ) , (2) 

where x = E^/En, (Jq = A 2 / 3 ijq (with <7o ^ 40 mb) provides a scaling of the 
cross-section with the atomic mass number (26), and F n (x, En) is a fragmentation 
function. Taking into account that inelastic hadronic collisions lead to roughly 
equal numbers of 7r°, ir + , and ir~ mesons, as well as the mild variation with 
pseudorapidity, the parametrization (28) 



1 + rxP(l -x?)J \l-xP l + rxP(l-xP) 



m 



\ 1/2 



1 - , (3) 

xE N ) y ' 

is found (29) to be consistent at la level with data collected at the Tevatron 
by the CDF detector (30); B w = a + 0.25, (3 = 0.98/^, r = 2.6/ y/a, a = 
3.67 + 0.83L + 0.075L 2 , and L = \n(E N /Te\). Substitution of Eq. (2) into 
Eq. (1) leads to 

where 



Q„(E n ) ~ Z AlT (a) Qi H (E w ) (4) 



Qf 1 (E N )=4cn H d ^(E N ) (5) 

and the spectrum-weighted moment of the inclusive cross-section or so-called 
Z-factor is given by 

l 

Z A n(a) = J ' X ' 1 F(x,E N )dx. (6) 
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On average, the photons carry one-half of the energy of the pion, and thus the 
7-ray emissivity is 

gmax 



Br n (s 7 ) 

~ Z n ,(a)Q^(E,), (7) 



where E™ m (E 7 ) = E 1 + m 2 r /(4£' 7 ) and 2' 7r7 (a) = 2/a. Assuming proton dom 



inance and emission spectrum oc E N 2 , the gamma ray emissivity above 1 TeV 



becomes 



Q 7 (> 1 TeV) = c (E^/En) a n H n CR (> 1 TeV) 
-29 / n H \ photons 



^ 10- 29 ) , (8) 

Vcm-^/ cm d s 

where (E^/En) ~ 0.2 is the average energy of the secondary pions relative to 
the cosmic ray protons. For different spectral indices, the quantity {E n /E^) is 
generalized to the spectrum-weighted moments for pion production by nucleons 
given in Eq. (6). 

The density of protons from a supernova converting a total kinetic energy W 
of 10 50 erg to proton acceleration is approximately given by p/W, where we will 
assume that the density in the remnant is not very different from the ambient 
energy density pe ~ 10 _12 ergcm -3 of galactic cosmic rays. This approximation 
is valid for young remnants in their Sedov phase (31). The total luminosity in 
gamma rays is then 

CJ> 1 TeV) = Q 7 — ~ 10 33 photons s _1 . (9) 
Pe 

The expected rate of TeV photons from a supernova at a distance d ~ 1 kpc is 



^ -(> 1 TeV) - ^ 



d(ln£ 7 ) v ' 4ird 2 

_ n ( photons \ / W 



~ 10" 



cm 2 s / \10 50 erg 



Neglecting absorption effects, the all flavor neutrino and photon fluxes are roughly 
equal, and so we anticipate an event rate of 3 detected neutrinos per decade of 
energy per km 2 yr, a result readily obtained from the relation 

dj\f u . m-ii ( neutrinos \ / area \ / time\ 



d(lnE u ] 



m , T . n / neutrinos \ / area \ / time \ „ . . 

< >1TeV) = 10 (^rj(rkrf)(T^J p -- 



where the last factor indicates the average probability for neutrino detection; for 
the TeV energy considered here P u ~ 10~ 6 (1). 

This estimate may be somewhat optimistic because we assumed that the neu- 
trino energy reaches 100 TeV with an E~ 2 spectrum. 1 On the other hand, if the 



Such a spectral index is typical of Fermi acceleration mechanisms (32). 
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galactic cosmic ray spectrum extends at least up to the knee (E ~ 3000 TeV), 
then some of the sources must produce ~ 100 TeV secondaries. 

Photo-disintegration of high-energy nuclei, followed by prompt photo-emission 
from the excited daughter nuclei also produces TeV gamma rays and antineu- 
trinos. In this chain reaction, the nuclei act in analogy to Einstein's relativis- 
tic moving mirror to "double-boost" eV starlight to TeV energies for a Lorentz 
boost factor > 10 6 (33). The TeV antineutrino counterpart is produced through 
/3-decay of the emitted neutrons (34). Effective nucleus photo-disintegration is 
expected in stellar associations harboring a large population of massive and early 
type stars that can provide the required Lyman-a emission. Directional beams 
of galactic antineutrinos and gamma rays may thus emerge as fluctuations over 
the diffuse intensity. 

A plethora of explanations have been proposed to address the production mech- 
anism of cosmic rays beyond the dip (35). In the absence of a single model which 
is consistent with all data, the origin of these particles remains a mystery. Clues 
to solve the mystery are not immediately forthcoming from the data, particularly 
since various experiments report mutually inconsistent results. In recent years, 
a somewhat confused picture regarding the nature of these cosmic rays has been 
emerging. The HiRes data have been interpreted as a change in cosmic ray com- 
position, from heavy nuclei to protons, at E ~ 4 x 10 8 GeV (36). This is an order 
of magnitude lower in energy than the previous crossover deduced from the Fly's 
Eye data (37). 

The end-point of the galactic flux is expected to be dominated by iron, as the 
large charge Ze of heavy nuclei reduces their Larmor radius (containment scales 
linearly with Z) and facilitates their acceleration to highest energy (again scaling 
linearly with Z) . The dominance of nuclei in the high energy region of the galactic 
flux carries the implication that any changeover to protons represents the onset 
of dominance by an extra- galactic component. The inference from this HiRes 
data is therefore that the extra-galactic flux is beginning to dominate the galac- 
tic flux already at E ~ 4 x 10 8 GeV. Such a transition would appear to require 
considerable fine-tuning of the shape and normalization of the two spectra. It has 
been argued however that the spectral shape required for the extra-galactic com- 
ponent can develop naturally during the propagation of extra-galactic protons 
in the cosmic microwave background (CMB) over cosmological distances (38). 
At E ~ 3 x 10 8 GeV the energy losses due to e + e~ pair production and cosmic 
expansion are roughly equal and thus produce a steepening of an initially feature- 
less power-law injection spectrum. Significantly, the onset of the extra-galactic 
component would be well below Eqzk ~ 5 x 10 10 GeV, the threshold energy for 
resonant P7cmb — ► A + — > Ntt energy-loss on the CMB, and so samples sources 
even at large redshift. 

It is helpful to envision the cosmic ray engines as machines where protons 
are accelerated and (possibly) permanently confined by the magnetic fields of 
the acceleration region. The production of neutrons and pions and subsequent 
decay produces neutrinos, gamma rays, and cosmic rays. If the neutrino-emitting 
source also produces high energy cosmic rays, then pion production must be the 
principal agent for the high energy cutoff on the proton spectrum. Conversely, 
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since the protons must undergo sufficient acceleration, inelastic pion production 
needs to be small below the cutoff energy; consequently, the plasma must be 
optically thin. Since the interaction time for protons is greatly increased over 
that of neutrons because of magnetic confinement, the neutrons escape before 
interacting, and on decay give rise to the observed cosmic ray flux. The foregoing 
can be summarized as three conditions on the characteristic nucleon interaction 
time scale Ti nt ; the neutron decay lifetime r ra ; the characteristic cycle time of 
confinement r cyc i e ; and the total proton confinement time r con f: (i) T m t 3> T cyc \ e ; 
(ii) T n > T cyc i c ; (iii) Tint <C T con f. The first condition ensures that the protons 
attain sufficient energy. Conditions (i) and (ii) allow the neutrons to escape the 
source before decaying. Condition (iii) permits sufficient interaction to produce 
neutrons and neutrinos. These three conditions together define an optically thin 
source (39). A desirable property to reproduce the almost structureless energy 
spectrum is that a single type of source will produce cosmic rays with a smooth 
spectrum across a wide range of energy. 

The cosmic ray flux above the ankle is often summarized as "one 3 x 10 10 GeV 
particle per kilometer square per year per steradian." This can be translated into 
an energy flux (41) 

3xl0 10 GeV 
^ JCR} (l io cm 2 ) (3 x 10 7 s) sr 

= 10~ 7 GeV cm -2 s -1 sr -1 . (12) 

From this we can derive the energy density pe in ultra-high energy cosmic rays 
using flux = velocity x density, or 

4vr J dE{EJ CR } = cp E . (13) 

This leads to 



R 

^max _ 

4tt f 10~ 7 Iri GeV iq TeV , „ 

PE = — / -E^dE ~3xl(T 19 ^, 14 

c J E cur cm 6 

F, ■ 

taking the extreme energies of the accelerator(s) to be E maiX /E mm ~ 10 3 . The 
power required for a population of sources to generate this energy density over 
the Hubble time of 10 10 yr is 



E 2 dN CR 



dE 



.[10 19 ,10 21 ] 



ln(10 21 /10 19 ) 



10 44 ergMpc _ V , (15) 



where an energy spectrum oc E~ 2 has been assumed (40). 

This works out to about 3 x 10 39 erg s -1 per galaxy, 3 x 10 42 erg s -1 per 
cluster of galaxies, ~ 2 x 10 44 erg s _1 per active galaxy, or ~ 2 x 10 52 erg per 
cosmological gamma ray burst (41). The coincidence between these numbers 
and the observed output in electromagnetic energy of these sources explains why 
they have emerged as the leading candidates for the cosmic ray accelerators. The 
coincidence is consistent with the relationship between cosmic rays and photons 
built into the previously discussed "optically thin" source. 
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The energy density of neutrinos produced through pj interactions of these 
protons can be directly tied to the injection rate of cosmic rays 

2 dN v 3 2 ^cR . , 

E »dK~8 e * tRE ^> (16) 

where tn is the Hubble time and e n is the fraction of the energy which is injected 
in protons lost into photo-pion interactions. (The factor of 3/8 comes from the 
fact that, close to threshold, roughly half the pions produced are neutral, thus not 
generating neutrinos, and one quarter of the energy of charged pion decays goes 
to electrons rather than neutrinos.) 2 The "Waxman-Bahcall bound" is defined 
by the condition e % = 1 



where the parameter £z accounts for the effects of source evolution with redshift, 
and is expected to be ~ 3 (42). For interactions with the ambient gas (i.e., pp 
rather than jry collisions), the average fraction of the total pion energy carried 
by charged pions is about 2/3, compared to 1/2 in the photo-pion channel. In 
this case, the upper bound given in Eq. (17) is enhanced by 33% (43). 

The actual value of the neutrino flux depends on what fraction of the proton 
energy is converted to charged pions (which then decay to neutrinos). To quantify 
this, we follow Waxman-Bahcall and define e n as the ratio of charged pion energy 
to the emerging nucleon energy at the source. For resonant photoproduction, the 
inelasticity is kinematically determined by requiring equal boosts for the decay 
products of the A + (44), giving e w = E n +/E n « 0.28, where E n + and E n are the 
emerging charged pion and neutron energies, respectively. For pp — > NN + pions, 
where N indicates a final state nucleon, the inelasticity is ~ 0.6 (45). This then 
implies that the energy carried away by charged pions is about equal to the 
emerging nucleon energy, yielding (with our definition) e n ~ 1. 

At production, if all muons decay, the neutrino flux consists of equal fractions of 
Ve, Vn and v„- Originally, the Waxman-Bahcall bound was presented for the sum 
of v„ and (neglecting v e ), motivated by the fact that only muon neutrinos 
are detectable as track events in neutrino telescopes. Since oscillations in the 
neutrino sector mix the different species, we chose instead to discuss the sum 
of all neutrino flavors. When the effects of oscillations are accounted for, nearly 
equal numbers of the three neutrino flavors are expected at Earth (46). (For a 
detailed calculation of the flavor ratios see e.g., (47).) 

If the injected cosmic rays include nuclei heavier than protons, then the neu- 
trino flux expected from the cosmic ray sources may be modified. Nuclei under- 
going acceleration can produce pions, just as protons do, through interactions 
with the ambient gas, so the Waxman-Bahcall argument would be unchanged in 

2 The average neutrino energy from the direct pion decay is (E V ^) 7T — (1 — r) E v /2 ~ 0.22 E n 
and that of the muon is {E fl } 7 ' = (1 + r) E^/2 ~ 0.78 E n , where r is the ratio of muon to the 
pion mass squared. In muon decay, since the has about 1 /3 of its parent energy, the average 
muon neutrino energy is (-E^ fl ) M = (1 + r)_E^/6 = 0.26 E n . 



[E 2 M W b ~ (3/8) £z e 7r £h -7— E 2 




w 2.3 x 10~ 8 e w i z GeV cm" 2 s" 1 sr 



-1 



(17) 
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this case. However, if interactions with radiation fields dominate over interac- 
tions with matter, the neutrino flux would be suppressed if the cosmic rays are 
heavy nuclei. This is because the photo-disintegration of nuclei dominates over 
pion production at all but the very highest energies. Defining k as the fraction 
of nuclei heavier than protons in the observed cosmic ray spectrum, the resulting 
neutrino flux is then given by 



Unfortunately, current observations do not allow a conclusive determination of 
k. For 1.6 < a < 2.1, the data can be well reproduced if the emitted cosmic 
rays consist entirely of nuclei with masses in the intermediate (carbon, nitrogen, 
or oxygen) to heavy (iron, silicon) range, but a mixture of protons (97%) and 
heavier species (3%) is also acceptable (16,48). 

The diffuse neutrino flux has an additional component originating in the en- 
ergy losses of ultra-high energy cosmic rays en route to Earth. Ultrahigh energy 
protons above the "GZK cutoff" interact with the cosmic microwave and infrared 
backgrounds as they propagate over cosmological distances. These interactions 
generate pions and neutrons, which decay to produce neutrinos (14). The accu- 
mulation of these neutrinos over cosmological time is known as the cosmogenic 
neutrino flux. Ultra-high energy nuclei also interact with the cosmic microwave 
and infrared backgrounds, undergoing photodisintegration. The disassociated 
nucleons then interact with the cosmic microwave and infrared backgrounds to 
produce cosmogenic neutrinos (49). In the limit that the cosmic backgrounds 
are opaque to cosmic ray nuclei, full disintegration occurs and the resulting cos- 
mogenic neutrino spectrum is not dramatically different from that predicted in 
the all-proton case (assuming the cosmic ray spectrum extends to high enough 
energies to produce protons above the GZK cutoff). In contrast, if a significant 
fraction of cosmic ray nuclei remain intact, the resulting flux of cosmogenic neu- 
trinos can be considerably suppressed. Of course, the predicted neutrino flux also 
depends on the spectrum, the source evolution, and the nature of the emitted 
cosmic rays (14-16). Representative spectra of cosmogenic neutrinos are shown 
in Fig. 1, for pure Fe and bi-modal Fe + p mass compositions consistent with 
data (16). The estimates shown in Fig. 1 are framed in the context of the usual 
concordance cosmology (50) of a flat universe dominated by a cosmological con- 
stant with 17a ~ 0.7, the rest being cold dark matter (CDM) with f2 m ~ 0.3. 
(A rigorous compilation of the ACDM parameters is provided in Sec. 4.5.) The 
Hubble parameter is given by H 2 (z) = Hq (S7 m (l + z) 3 + Qa), normalized to its 
value today of 70 kms -1 Mpc -1 . The time-dependence of the red-shift can be ex- 
pressed via dz = —dt (1 + z)H. The cosmological evolution of the source density 
per co-moving volume is parameterized as 



where TC(z) is the cosmological evolution of the cosmic ray sources, which is taken 
here to follow the luminosity density evolution of QSOs: T~Iqso{z) = (1 + -z) 3 , for 
z < 1.9, Hqso{z) = (1 + 1-9) 3 for 1.9 < z < 2.7, and H{z) = (1 + 1.9) 3 exp{(2.7 - 
z)/2.7}, for z > 2.7 (15). For such a cosmological evolution, which is similar to 
that describing the star formation rate, £z ~ 3 (42). 



El J u « (1 - k) [£*$„]wb • 



(18) 



C l (z,E)=H{z)£ i {0,E), 



(19) 
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In summary, the case for doing high energy neutrino astronomy is compelling; 
the challenge is to deliver the technology to build a neutrino detector with the 
largest possible effective area and the best possible angular and energy resolution 
sensitive to these cosmic neutrino fluxes. We discuss this next. 

3 NEUTRINO TELESCOPES 
3.1 Detection principle 

It took about 35 years between Markov's conception of the Neutrino Telescope 
(NT) detection principle (51) to the operation of the first full detector. Though 
NTs use the well established technique of photomultipliers (PMTs), they are 
challenging instruments to build. The difficulties for construction arise from (i) 
the necessity of instrumenting extremely large regions, mainly because of the 
smallness of the neutrino-nucleon cross section and the steeply decreasing energy 
spectrum of neutrino fluxes of atmospheric and astrophysical origin; (ii) the re- 
quired darkness and transparency of the medium to detect the faint Cherenkov 
light produced by ultra-relativistic charged particles; and (Hi) the need for filter- 
ing neutrino events out of the more abundant atmospheric muons. Detectors are 
located under 1—4 km of water or ice, where, for E v > 100 GeV, the downward- 
going background of atmospheric muons is about 5 — 6 orders of magnitude larger 
than the atmospheric neutrino flux coming from all directions (see Figs. 2 and 
3). 

NTs are tridimensional matrices of photodetectors made of strings of optical 
modules (OMs) that house large photocathode PMTs and protect them from 
both the water column pressure (when they are installed in the sea or in lakes) 
and the pressure during ice refreezing (for Antarctic detectors). In some cases 
OMs also comprise front-end electronics. The NT energy threshold depends on 
the vertical spacing between PMTs and the horizontal distance of strings; typi- 
cal values range between 50 GeV to 1 TeV. The threshold is not sharp, but the 
detection efficiency increases with energy. Neutrinos interact with the medium 
nuclei in and around the detector through charged (z-'iepton + X — ► lepton + X) 
and neutral (fi cp ton + X — > lepton + X) current interactions. Charged current 
deep inelastic collisions above 10 GeV dominate over all other processes with the 
exception of v e e interactions, because of the intermediate-boson resonance formed 
in the neighborhood of ££ es = M^/2m e « 6.3 x 10 6 GeV, generally referred to 
as the Glashow resonance (62). The resonant reactions v e e~ — ► W~ — > Vafi~ an d 
v e e~ — > W~ — > hadrons may offer a detectable signal if the v e flux originates 
in pp collisions (43). 3 The neutrino-nucleon cross section rises almost linearly 
with energy (63). This implies that (i) the interaction probability increases with 
energy, indicating that NTs are optimal detectors at high energies, and (ii) neu- 
trinos begin to be absorbed in the Earth. The shadowing effect of the Earth not 
only depends on the neutrino energy but also on the incident zenith angle (64). 

3 In pp collisions the nearly isotopically neutral mix of pions will create on decay a neutrino 
population in the ratio N Vfl — N^^ — 2N Ve — 2N„ C . In contrast, photopion interactions leave 
the isotopically asymmetric process p7 — » A + — » ir + n as the dominant source of neutrinos; then, 
at production, N Vfi = N Vfl = N Ve > N„ e . 
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At about 10 5 GeV, the interaction length of muon and electron neutrinos is 
about the Earth's diameter, and hence about 80% (40%) of the neutrinos with 
cos# = — 1 (—0.7) are absorbed, where 6 is the zenith angle. Hence, this indi- 
cates that at very high energy, where it is expected that the astrophysical signal 
emerges from the atmospheric muon and neutrino backgrounds, detectors should 
have a large acceptance. There is a completely different behavior for tau neutri- 
nos. In propagating through the Earth, they interact producing tau-leptons that 
decay back into a v T of lower energy in what is referred to as a "regeneration 
effect" (65). The two leptonic channels, with a branching fraction of about 17%, 
produce also v e and u^. The relevance of this phenomenon in terms of event rates 
at detectors depends on the neutrino spectrum. 

NTs detect ultra-relativistic (j3 ~ 1) secondary particles produced in the inter- 
action of the incoming neutrino with an atomic nucleus. These particles, which 
typically travel faster than the speed of light in the detector medium, emit di- 
rectional Cherenkov light at an angle 9c = cos -1 [l//3n] ~ 41°, where n is the 
refraction index of the medium. Between 300-600 nm, the region where PMTs are 
mostly sensitive, about 3.2 x 10 4 Cherenkov photons are emitted by a bare muon 
track in the ice per meter. Once the OM acceptance and PMT quantum efficiency 
are accounted for, on average only 6% of these photons contribute to the signal. 
The muon neutrino event topology is very different from that of electron and 
tau neutrinos. Secondary muons travel long distances and can be reconstructed 
with an angular resolution better than a degree at high energy. Namely, the 
kinematic angle between the muon and the neutrino in the deep inelastic regime 
(> 10 GeV) decreases with neutrino energy as y / '(O^) oc \J m p /E u (26), with m p 

the proton mass. For E u > 10 TeV, the angle becomes negligible with respect 
to the intrinsic detector resolution. NTs measure the leading edge time of the 
PMT signal and in some cases the fully digitized waveform and the integrated 
charge. This information is used to reconstruct tracks for muons that travel long 
distances across the detector and determine their energy as well. Charged cur- 
rent (CC) interactions of electron and tau neutrinos and all flavor neutral current 
(NC) interactions produce cascade-like events. The resulting hadronic and elec- 
tromagnetic showers are small with respect to the scale of the separation between 
OMs, and so they are almost a point- like source of light that propagates with 
some anisotropy that remembers the direction of the incoming neutrino. This 
prevents the achievement of an angular resolution for neutrinos comparable to 
the muon channel. 

Muons lose energy according to 

where a = 2.0 xlO" 6 TeV cm /g is the continuous ionization term. The stochastic 
energy loss term bE, with b = 4.2 x 10~ 6 cm 2 /g (50), dominates at energies 
> 1 TeV in water or ice. The distance a muon travels before its energy drops 
below some energy threshold Efi, called the muon range, is then given by 



A, = ^ln 



a + bE^ 
a + bE^ 



(21) 
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For muons with > 1 TeV propagating in water or ice, stochastic losses, such 
as bremsstrahlung, pair production, and photonuclear processes (bE^ term) dom- 
inate over the ionization constant a term. In the first kilometer, a muon of about 
100 TeV typically loses energy in 2 — 3 showers carrying more than 10% of its 
initial energy. Near the end of its range the muon becomes a minimum ionizing 
particle emitting light that creates single photoelectron signals at a distance of 
just over 10 m from the track. Above 1 TeV, the muon energy can be recon- 
structed with an energy resolution in log 10 £^ of about 0.3. Of course, above 
a certain energy saturation effects of the PMT electronics and the size of the 
detector limit the energy reconstruction ability. In fact, because many of the 
events are not fully contained at very high energy, part of the Cherenkov cone 
can miss the PMTs so that part of the charge, that is proportional to the muon 
energy, is lost. The parent neutrino energy can be inferred from simulation of the 
v{y) + N — > fi~(fi + ) interactions. For E v ^ > 10 5 GeV, the average value of the 
energy taken by the lepton is ~ 0.75 (63). In most cases, the well reconstructed 
cascade events develop inside the detector and hence the neutrino energy reso- 
lution is largely improved, up to 0.1 - 0.2 in log 10 ^. Tau neutrinos produce 
similar events to electron neutrinos up to about 10 6 GeV, because it is not yet 
demonstrated that NTs have the ability to discriminate between electromagnetic 
and hadronic showers. However, for B i; ~2x 10 6 GeV, once again tau neutrinos 
deserve a special attention: for such an energy, the tau-lepton range is about 
100 m and consequently both the v T N interaction and the shower (hadronic or 
electromagnetic) from the tau decay can be separated (66). These are called 
double-bang events and are in principle a background-free topology. Neverthe- 
less, the expected event rates are limited by the small energy window in which 
these events are detectable. As a matter of fact, this topology is limited in en- 
ergy to the range of a few PeV for the tau track to be long enough, up to about 
100 PeV to contain the two cascades in the instrumented volume (the r range 
is 1 km at about 200 PeV). The three different neutrino topologies are shown in 
Fig. 4. 

Atmospheric muons and neutrinos constitute a troublesome background for 
searches of neutrinos produced in astrophysical sources. Moreover, the predic- 
tions for such a background are quite uncertain at energies > 1 TeV. Therefore, 
background discrimination at NTs requires different techniques depending on the 
analysis. However, all of these analyses use the fact that neutrinos can cross the 
entire Earth, though with a shadowing effect that depends on energy and nadir 
angle. By selecting upgoing muons induced by neutrinos, the residual background 
of atmospheric muons is largely reduced to those events that are downgoing but 
are mis-reconstructed as upgoing. For cascade-events background discrimination 
is instead accomplished through containment cuts, since neutrinos have their 
vertex in the instrumented region while external tracks are produced mostly by 
atmospheric muons. 

For point source searches, the directional information is the crucial cipher. 
Being undeflected by intergalactic and galactic magnetic fields, neutrinos preserve 
directionality from their own sources and tend to cluster around it while the 
background distribution is random (see Sec. 4.1). The golden channel to point 
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back to sources is the muon one. This is because the angular resolution is better 
than 1°, whereas cascade like events have poor resolution of about 10° — 30°, 
depending on the medium. 

For diffuse flux searches, the energy measurement is determinant. This is 
because the signal spectrum is expected to be harder than the neutrino back- 
ground (95) from atmospheric cascades initiated by cosmic rays. When protons 
and nuclei enter the atmosphere, they collide with the air molecules and produce 
all kind of secondary particles, which in turn interact or decay or propagate to 
the ground, depending on their intrinsic properties and energies. Above 1 GeV, 
the most abundant particles at sea level are neutrinos and muons (50). The prop- 
erties of these particles, which constitute the main background for astrophysical 
searches, are described in Sec. 4.2. 

3.2 Detectors in operation and decommissioned 

The DUMAND and the Baikal collaborations pioneered underwater NT technolo- 
gies in the deep ocean near Hawaii (67) and in the Siberian Lake Baikal (68), re- 
spectively. In November 1987, the DUMAND Collaboration measured the muon 
vertical intensity at depths ranging between 2 — 4 km, with a prototype line at a 
depth of 4.8 km about 30 km off-shore of Hawaii (55). The project was canceled in 
1995. In 1993 a first configuration of 36 PMTs on 3 strings (NT-36) was installed 
in Lake Baikal in Siberia at the shallow depth of 1.1 km and 3.6 km off-shore (54). 
The experiment was later upgraded to larger configurations (NT-72-96-144). In 
700 days of effective livetime, 320 million muon events were detected, and the 
first atmospheric neutrino measured by a neutrino telescope was selected in this 
sample and reported in 1996. The configuration of the experiment (NT-200) with 
192 OMs was put into operation in April 1998. It consists of an umbrella-like 
structure of 8 (72 m long) strings, with up-looking and down-looking OMs (con- 
taining 37-cm diameter PMTs developed for the project, called QUASAR-370). 
Three external strings at 100 m from the center of NT200, each with 12 pairs of 
OMs, were added in April 2005 to increase the cascade sensitivity at very high en- 
ergies (NT200+). The Baikal Collaboration is working towards the construction 
of a Gigaton-volume detector with about 2100-2500 OMs (arranged over 95-100 
strings separated by distances of the order of 80-120 m) grouped in clusters to 
form independent sub-arrays (69). A prototype string has been taking data since 
April 2008. 

The NESTOR Collaboration began surveys of an area close to the Pelopon- 
nese coast at a mean depth of about 4 km in the '90s. They measured the 
muon flux with a hexagonal prototype floor with 12 up-looking and down-looking 
PMTs (56), and studied by simulation the possibility to deploy a telescope made 
of towers with 12 of these hexagonal floors (32 m in diameter) vertically spaced 
by 20 — 30 m. The Collaboration is operating a deep sea station off-shore of 
Pylos at a depth exceeding 3.5 km connected to the on shore laboratory by an 
electro-optical cable. 

In 1987 Halzen came out with the idea of "ice-fishing" neutrinos. The Antarctic 
Muon And Neutrino Detector Array (AMANDA) was located below the surface 
of the Antarctic ice sheet at the geographic South Pole (70). During 1993 and 
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1994, in an exploratory phase, the four-string AMANDA — A array was deployed 
and instrumented with 80 PMTs spaced at 10 m intervals from 810 to 1000 m. 
(The scattering length at that depth is too short to allow degree level reconstruc- 
tion.) A deeper array of 10 strings, referred to as AMANDA-B10, was deployed 
during the austral summers between 1995 and 1997, to depths between 1500 and 
2000 m. The instrumented volume of AMANDA-B10 formed a cylinder with 
diameter 120 m, overlooked by 302 optical modules (Hamamatsu R5912-2 20 cm 
diameter) (71). During December 1997 and January 2000, the detector was ex- 
panded by adding nine outer strings of PMTs. The composite AMANDA-II array 
of 19 strings and 677 OMs forming two concentric cylinders with larger diameter 
of 200 m became operative in 2000 and continued up to 2009. The OMs were 
connected to cables supplying the 8-inch PMT power and transmitting analog 
signals. Some of the strings used fiber optic cables reducing cross-talk effects, 
but most used coaxial twisted-pair cables. One of the strings pioneered the digi- 
tization technique, with data acquisition for waveform digitization in each module 
communicating to the surface, designed as a prototype for the IceCube Digital 
OMs (DOMs) and data transmission. 

The IceCube neutrino telescope, currently being deployed near the Amundsen- 
Scott station, comprises a cubic-kilometer of ultra-clear ice about a mile below 
the South Pole surface, instrumented with long strings of sensitive photon de- 
tectors which record light produced when neutrinos interact in the ice. The 
deep array is complemented by IceTop, a surface air shower detector consisting 
of a set of 160 frozen water tanks. With its surface array, IceCube becomes a 
3-dimensional air shower array for studies of cosmic rays up to EeV. It is also 
useful for calibration of the neutrino telescope, and it serves as a partial veto for 
the downward events at high energy. IceCube strings are placed 124 m apart on 
a hexagonal grid. Each string is instrumented (between 1.5 to 2.5 km below the 
surface) with 60 DOMs vertically separated by about 17 m. The DOMs enclose 
10-inch Hamamatsu PMTs with quantum efficiency ~ 20% at around 400 nm. 
Six additional strings with new high efficiency (~ 25 — 30% at 400 nm) PMTs 
to be deployed in between 7 standard IceCube strings will constitute the Deep 
Core. It will improve the sensitivity of the detector in the region < 1 TeV, im- 
portant for dark matter searches and neutrino oscillations. The outer IceCube 
strings around Deep Core and the upper layers of DOMs will be used as veto- 
shield against the background of atmospheric muons allowing the identification 
downgoing neutrino induced muons when the neutrino vertex is in the inner de- 
tector. Construction of the observatory will be completed in the austral summer 
2010-11. The drilling of strings has been so successful that up to about 20 strings 
can be installed in a season (from mid- November to mid-February, when airplane 
flights are possible to and from the South Pole). The basic element of IceCube, 
the DOM, is made of a glass spherical pressure vessel of 33 cm diameter, con- 
taining a 25 cm diameter PMT, a high voltage supply, light emitting diodes for 
energy and time calibrations, and an electronic motherboard. The motherboard 
digitizes and time-stamps PMT signals (72). The performance of the first string 
(deployed in the 2005-6 season), its calibration and test methods, as well as the 
data acquisition and installation of the detector are described in (73). The layout 
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of IceCube is shown in Fig. 5. 

The completion of the undersea neutrino observatory ANTARES, located at 
about 2.5 km below the sea surface off-shore Toulon in South France, took several 
years. Since 1996, the Collaboration had sea campaigns comprising more than 
60 line deployments for tests of water transmission, light optical background, 
bio-deposition and sedimentation (74). The permanent electro-optical cable, of 
about 40 km, transmitting data and power between the shore station in La Seyne 
sur Mer and the junction box off shore has been operating since November 2002. 
Between 2003 and 2005 various prototype lines and a mini-instrumentation line 
for environmental parameter measurements were deployed (75). In 2006, the 
first 2 lines of the detector were deployed and 8 additional lines were installed 
in 2007. On May 2008 ANTARES was completed with a total of 12 lines. Lines 
were connected initially by a submarine and then by an unmanned Remote Op- 
erated Vehicle. The ANTARES observatory comprises the 12 mooring lines, a 
line specifically dedicated to marine environmental monitoring, a seismometer, 
and a biocamera for bioluminescence studies. The lines are anchored to sea bed 
at 2475 m depth and held vertical by buoys. Buoys are freely floating so each 
line moves under the effect of the sea current, with movements of a few meters 
for typical values of 5 cm/s. An acoustic positioning system, made of transpon- 
ders and receivers, gives a real time measurement of the position of the OMs 
with a precision better than 10 cm, typically every 2 minutes. Tiltmeters and 
compasses provide their orientation. Seventy-five OMs along each lines between 
about 2400 and 2000 m are grouped in triplets on storeys. The 3 PMTs (Hama- 
matsu 10" (76)) look downward at 45° from the vertical to prevent transparency 
loss due to sedimentation and optimize detection for upgoing tracks. PMTs are 
housed inside pressure resistant glass spheres made of two halves closed by ap- 
plying an under-pressure of 200-300 mbar. The set-up, including the PMT, the 
glass sphere, the silicon gel for optical coupling between the glass and the photo- 
cathode, and the mu-metal cage for shielding the Earth magnetic field is referred 
to as the OM (77). Storeys also include titanium containers housing the front- 
end electronics, with a pair of ASIC chips per PMT used for signal processing 
and digitization. This provides the time stamp and the amplitude of the PMT 
signal. Each of the OMs contains a pulsed LED for calibration of the relative 
variations of PMT transit time and a system of LED and laser Optical Beacons 
allows the relative time calibration of different OMs. An internal clock system, 
which is synchronized by GPS to the Universal Time with a precision of ~ 100 ns, 
distributes the 20 MHz clock signal from the shore. Time calibrations allow a 
precision at the level of 0.5 ns ensuring the capability of achieving an angular res- 
olution at the level of 0.3° muons above 10 TeV for point source searches (78). 
All data above a threshold of about 1/3 of a photoelectron pulse is sent to shore 
for further online filter. This requires coincidences between PMTs on the same 
storey rejecting hits not compatible to the propagation of light in water between 
hit PMTs (79). The ANTARES detector is shown in Fig. 6. 

Ice and water are proving to be suitable transmission media with different pros 
and cons. Photon propagation properties are characterized by a much longer 
scattering length A scat t in sea water compared to ice. The effective scattering 
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length, denned in terms of the average scattering angle of the photons, 9, reads 
A e flf = A SC att/(l — (cos 9}). In sea water, A c g ~ 100 — 200 m is much larger than 
in ice due to a more forward peaked angular distribution. On the other hand the 
absorption length, that determines the distance photons can propagate before 
being absorbed, is much larger in ice (~ 100 — 200 m) than in water ~ 50 — 70 m. 

Sea water is a less quiet environment than ice. Additionally, bioluminescence 
and /3-decay of 40 K can produce a baseline of optical background in the sea of 
the order of 30 — 100 kHz and MHz bursts that cause a random dead time for 
OMs. No current or bioluminescence are present in ice. Once water developed 
in the drilling process is frozen, the risk of short circuits due to water leakages in 
the OMs is largely reduced, and PMTs can be switched on. On the other hand, 
once ice freezes, the components are not recoverable, while entire instrumented 
lines can be recovered from water for repairs. The average background rate in 
IceCube PMTs is ~ 280 Hz, mainly dependent on the radioactivity from the OM 
and PMT materials. The lower optical noise in ice makes IceCube a very sensitive 
detector for SN collapse searches. 4 

The neutrino effective area, the sensitive area 'seen' by neutrinos producing 
detectable muons when entering the Earth, is a useful parameter to determine 
event rates and the performance of a detector for different analyses. The expected 
event rate 



of a model predicting a neutrino spectrum J v depends on the neutrino-nucleon 
cross section a(E u ), the shadowing of the Earth PEa,rth(E u ,9 u ), the track recon- 
struction quality cuts, and the selection criteria for background rejection. The 
shadowing factor depends on column depth, that is the integral of the Earth den- 
sity profile p(9 u , I) for a given direction of the neutrino 9 V over the distance in the 
Earth L P Ea rth{K, V ) = e -NAXv(E„)xf e p(0„,e)M^ where Na ig the Avogadro num- 
ber. The neutrino effective area is a parameter calculable through simulations 
and formally given by 

Af(E v , Q u ) = V e s(E l/ , 9 V , <j> v ) x (pN A ) x <j{E v ) x P Earth (E u , 9 U ) (23) 

where (pNX) is the number of nucleons in the target material assumed isoscalar 
and V cS (E u ,9 u ,phi u ) = V gen (E u ,6 u ,(f) u ) x N SC \(E U , 9 U , 4> u ) / N gcn (E u , U , (f>„) is the 

4 About 99% of the binding energy of a Type II SN is expected to be released in neutrinos. 
These SNe originate in the gravitational collapse of massive giant stars, heavier than 8 solar 
masses, into a neutron star. Current upper limits on the number of supernova explosions in 
our Galaxy, from scintillator and water Cherenkov detectors, are getting closer to the predicted 
rates of about 2 per century (85). In the early phase (~ 10 ms), called deleptonization or 
neutronization, electron neutrinos are emitted in the electron capture process [e~ +p — » n + v e ) 
and subsequently (~ 10 s) all neutrino flavors are produced in reactions such as e + + e~ — » v + v 
in the thermalization phase. The passage of a large flux of MeV-energy neutrinos during a 
period of seconds will be detected as an excess of single counting rates in all individual optical 
modules. Each OM will be triggered by positrons from the inverse-/? reaction: v e +p — > e + + n, 
where the visible positron energy is given by -E v is = E De — Q + m e = E De — 0.789 MeV (86). 
IceCube, with 5160 optical modules, can observe the 10 ms pulse as an increased rate over the 
average background rate (280 Hz) in the PMTs. For example, a SN at a distance of 7.5 kpc 
would produce a rate of about 10 6 Hz over this time, yielding a 34a effect. A SN in the Large 
Magellanic Cloud would produce a 5cr effect. 




(22) 
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effective volume of the interactions generated in a volume V gcn of the neutrinos 
that produce detectable muons. The effective volume depends on the energy of 
the neutrinos that determines how far muons can travel from the vertex (muon 
range) and accounts for the selection efficiency of an analysis (the efficiency factor 
N se \/Ng Cn ). The effective areas of the different NTs are compared in Fig. 7. 

The point spread function (PSF) is the percentage of muons that are recon- 
structed inside an angular half-cone width from the neutrino or source direction. 
It depends strongly on the medium properties and on the reconstruction algo- 
rithms. In the final configuration of IceCube, the angular resolution at energies 
between 10 — 100 TeV will be such that about 50% of the reconstructed muons 
from the direction of a neutrino point source will be inside an angle of 0.5°. 
Sea water has longer effective scattering lengths than Antarctic ice hence the 
resolution in ANTARES or in a future cubic-kilometer detector in the Mediter- 
ranean (84) will be of the order of 0.2° — 0.3°. The PSF of IceCube is shown in 
Fig. 8. 

Cosmic ray experiments, like the Pierre Auger Observatory, provide a com- 
plementary technique for ultra-high energy neutrino detection by looking for 
deeply-developing, large zenith angle (> 60°) air showers (87). At these large 
angles, hadronic air showers have traversed the equivalent of several times the 
depth of the vertical atmosphere and their electromagnetic component has been 
extinguished far away from the detector. Only very high energy core-produced 
muons survive past 2 equivalent vertical atmospheres. Therefore, the shape of a 
hadronic (background in this case) shower front is very flat and very prompt in 
time. In contrast, a neutrino shower appears pretty much as a "normal" shower. 
It is therefore straightforward to distinguish neutrino induced events from back- 
ground hadronic showers. Moreover, very low v T fluxes could be detected very 
efficiently by Auger detectors by looking at the interaction in the Earth crust of 
quasi horizontal v T inducing a horizontal cascade at the detector (88). 

4 EXPERIMENTAL RESULTS AND SCIENTIFIC IMPACT 

The main goal of neutrino telescopes is detecting astrophysical neutrinos with 
E v > 100 GeV. As we noticed in Sec. 3, muons constitute the "golden sample" 
for neutrino astronomy since the achievable angular resolution is at sub-degree 
level. At energies > 10 TeV the limiting angle is mainly due to the transit time 
spread of the PMTs and to the scattering of light. 

While Imaging Atmospheric Cherenkov Telescopes (IACTs) can prove their 
point spread function on copious sources of gammas, such as the Crab Nebula, 
no astrophysical high energy neutrino source has proved to exist yet. Hence, the 
'standard candle' for neutrino telescopes must be an 'anti-source', such as the 
Moon shadow. The Moon disk, with a diameter of about 0.5°, blocks primary 
cosmic rays hence producing a deficit in the muon flux. IceCube, using 3 months 
of data of the 40 string configuration, begins to find evidence at 5a level of this 
deficit, which is a proof of the absolute timing and pointing capability of the 
detector (89). A detector with an ice platform on top (like IceCube) is able to 
test, for a limited solid angle, the angular resolution using events in coincidence 
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with the extensive air shower array IceTop. Sea detectors, such as ANTARES, 
can use instead a scintillator array carried by a boat, or the reconstruction of the 
shape of the close-by continental crust using horizontal muons (90). 

4.1 Steady and flaring point-source searches 

Searches of point-like sources have various advantages when compared to searches 
of diffuse fluxes, as in the former one can use the direction, the energy, and eventu- 
ally the time distributions to reject the background of atmospheric neutrinos. As 
a matter of fact, the signal is expected to distribute around the source according 
to the detector PSF, and to have a harder energy spectrum than the background. 
Flaring sources as well as GRBs are particularly interesting, because if the neutri- 
nos are emitted in coincidence with the gamma flare it is possible to dramatically 
reduce the background of neutrinos over a short time window. Since the back- 
ground of the experiment is fiat in right ascension (which is directly connected to 
the hour of the day), many 'equivalent experiments' with only background can be 
reproduced by scrambling the right ascension of events and keeping fixed other 
information, such as the declination and the energy proxy fixed. This has the 
advantage that data can be used for estimating the background rather than using 
simulation, so that resulting significances of the signal and the background are 
not affected by a possible imperfect understanding of the detector, or similarly 
the detector medium. Moreover, for searches of flares or GRBs, data outside the 
short time window can be used to reproduce 'equivalent experimental' samples 
where no signal is present. These techniques are similar to the off-on source tech- 
niques applied in photon astronomy, where the background outside the source 
pixel region is used to evaluate the significance of the observed source. In the 
many neutrino astronomy point-source searches (91), binned methods that would 
look for clusters of events on top of the background from specific directions or 
in the entire sky were applied. The search bin was optimized depending on the 
expected signal shape (typically assumed to be oc E~ 2 ) but to avoid missing 
signal events at the border of bins, grids or bins had to be moved, hence, penal- 
izing the search with many trial factors. Many experiments (92) use likelihood 
methods (93) that can be up to 40% more sensitive depending on the amount 
of information used. These methods use probability distribution functions that 
maximize the discrimination between the signal and background. The signal is 
characterized with respect to the background by the source directional feature 
and by the harder spectrum expected from astrophysical sources compared to 
atmospheric neutrinos. For the case of time dependent sources, such as periodi- 
cal binaries like micro-quasars or flares from blazars, periodicity assumptions or 
light curves from gamma, X-ray, and optical telescopes can be used. 

The IceCube's 22 string configuration collected 5114 events from the North- 
ern hemisphere in 276 days (80). No signal has been found and upper limits 
have been set for sources in the Northern hemisphere. For illustrative pur- 
poses we note that, for Geminga, the Crab Nebula, Markarian 421, M87, and 
Cygnus OB2, the upper limits of the 90% CL interval for a flux $ < <£go x 
lO-^TeV^cm-V^/TeV)- 2 are $ 90 = 9.67, 10.35, 14.35 , 7.91, 15.28, re- 
spectively. In addition, a hot spot was found with a post-trial probability of 
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1.34%. The excess, however, was not confirmed using the 6 month data of the 
40 string configuration (that by then encompassed in sensitivity the previous 
configuration by about a factor of 2) and so the hot-spot has been considered 
a background fluctuation (83). This 40 string analysis selects upgoing neutrino 
events from the Northern hemisphere at a rate of about 40 per day. In this re- 
gion of the sky the experiment is sensitive mainly to events in the 10 — 100 TeV 
region for a oc E~ 2 spectrum. Very recently, by using scrambled background 
of very high energy atmospheric muons, the IceCube Collaboration extended 
the field of view for point source searches to the Southern hemisphere. In this 
area of the sky, the sensitive region, for a spectrum oc E~ 2 , moves to higher 
energies (PeV-EeV) (94). As an illustration, the resulting 90%CL upper limits 
$ < $90 x 10~ 9 GeV" 1 cm" 2 s" 1 (E u /GeV)- 2 , using 22 strings of IceCube, 
are $ 90 = 558.2, 245.0, 27.1, for Centaurus A, Sagittarius A*, and 3C279 re- 
spectively. (The associated energy intervals that contain 90% of the signal are: 
AEv/GeV = 2 x 10 6 - 8 x 10 8 , 1 x 10 6 - 9 x 10 8 , 5 x 10 4 - 1 x 10 8 .) For the full 
data sample of 40 strings, the IceCube discovery reach (5a) for E 2 §go (averaged 
over right ascension, and in the declination region between 0° — 90°) is expected 
to be between 1 to 5 x 10~ n cm~ 2 s -1 . This is near to the exclusion flux for 
22 strings and about a factor of 3 higher than the exclusion limits for the same 
configuration. With the 40 strings result, flux limits for galactic sources are be- 
ginning to be in the region of interest: E 2 $ ~ 10~ 12 — 10 -11 TeV cm~ 2 s _1 (see 
Sec. 2). Actually, according to estimates based on gamma- ray observations it 
would be possible to achieve 5cr significance for Milagro Pevatrons, with ~ 10 yr 
of IceCube data collection (6). The 40 string sky-map is shown in Fig. 9. 

4.2 Atmospheric muons and neutrinos 

Atmospheric muons and neutrinos detected by NTs above some tens of GeVs are 
produced in cosmic ray induced showers and come mainly from pions and kaons. 
At lower energies (< 10 GeV) muon and electron neutrinos are also generated by 
muon decays. They were extensively studied in the past two decades because of 
the so called "atmospheric neutrino problem," a deficit of muon induced neutrino 
events that was explained by muon neutrino oscillations into tau neutrinos [for 
a review see, e.g., (95,96)]. At energies above 100 GeV, where NTs operate, 
atmospheric neutrino oscillation studies are not a primary target. Nonetheless, 
although DeepCore and ANTARES photocathode coverage are modest, they will 
collect a very large statistics of atmospheric neutrinos opening the possibility of 
observing neutrino oscillations. The main interests for NTs to measure atmo- 
spheric neutrinos are (i) calibration of the detector, demonstrating their ability 
to discriminate neutrino events from the huge background of atmospheric muons; 
(ii) understanding of hadronic interactions in the forward region, at energies not 
yet accessible by colliders; and (in) study of new physics scenarios. 

The dominant production channels of atmospheric and muon neutrinos de- 
tected in NTs are the decay chains of mesons created in hadronic cascades. The 
pion decay chain (ir + — ► — > e + v e v lA p^ L and the conjugate process) domi- 
nates the atmospheric neutrino production at ~ GeV energies. At energies of 
some GeVs, the muon decay length begins to be comparable to the atmospheric 
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altitude where most rrmons and neutrinos are produced (10 — 20 km). Hence, 
depending on the incoming zenith angle, muons may stop decaying. At ener- 
gies of about 100 GeV kaon decay becomes increasingly important with respect 
to pion decay. In this energy range, the atmospheric neutrino spectrum can be 
parametrized as (26) 



d JK _ A E , 

d\nE u tot " [l + B p E„(cos6)*/e n ' i + B k E u (cos6)* / e k 



A " + 0.635 x C " 



(24) 



where 0.635 is the branching ratio of decays; = 115 GeV and ex = 
850 GeV are the critical energies respectively for pion and kaon decay, a ~ 1.7 is 
the spectral index of the primary spectrum, Ai ot ] A v and B v are normalization 
constants; and B p and B k are constant terms with energy. The term (cos 6)* is 
equal to cos 6 for zenith angles such that cos 9 < 0.3, and for horizontal directions 
(cos 8 > 0.3) it has a larger value that accounts for the curvature of the Earth. 
The 1/(E U cos 6) dependence, is characteristic of atmospheric processes where 
there is competition between interaction and decay, since the decay probability 
increases with respect to the interaction one for longer horizontal paths in the 
atmosphere. The muon flux in this energy region has a similar dependence with 
the important difference on the energy fraction taken by muons in two body 
pion and kaon decays. We remind the reader that, on average, in pion decay 
the neutrino takes a quarter of the energy of the pion while the muon takes 
the rest. This implies that above ~ 1 TeV, muon production is dominated by 
pion decay and neutrino production by kaon decay (about 85% of the muons are 
produced via pion decay, whereas only 25% of neutrinos come from the decay 
of charged pions). The electron neutrino formula is more complicated, with two 
similar terms due to and Kl decays from which v e are produced directly, 
and a term due to muon decay into electron neutrinos determined by ir^, 
and Kl decays. In this energy window, the flavor ratios are w e : Wu : w T ~ 
1/20 : 19/20 : (97) and the energy spectra are a function of the zenith angle of 
the atmospheric cascades. As already observed, mesons in inclined showers spend 
more time in rarefied atmosphere where they are more likely to decay rather than 
interact. For this reason the spectra of highly inclined neutrinos are harder than 
those of almost vertical neutrinos and resemble their primary spectrum. Above 
about 10 5 GeV, kaons are also significantly attenuated before decaying and the 
"prompt" component, arising mainly from very short-live charmed mesons (D^, 
D°, D s and A c ) dominates the spectrum (18). Such a "prompt" neutrino flux is 
isotropic with flavor ratios 12/25 : 12/25 : 1/25. Since prompt neutrinos have a 
hard spectrum resembling the primary one (below the knee the differential cosmic 
ray spectrum goes as ~ E~ 2 - 7 ), they are a persisting background for astrophysical 
neutrinos that follow a power law ~ E~ 2 . 

Various Monte Carlo high precision calculations have been developed for study- 
ing the atmospheric neutrino oscillations (58,95,98). For energies below 10 GeV, 
high precision needs to account for the angle of secondaries with respect to pri- 
maries (99), the solar modulation effects, and the Earth geomagnetic field de- 
flections. At higher energies these effects become negligible and one dimensional 
calculations that neglect the angle of secondaries respect to primaries are suit- 
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able. Uncertainties in atmospheric neutrino fluxes are of the order of 15% in the 
region of ~ 10 GeV, rising to about 20-30% at ~ 100 GeV-1 TeV (100). These 
are mainly due to the rather poor knowledge of the primary cosmic ray flux and 
the hadronic interaction models. A precise estimate in the prompt region is not 
possible since charm production in the atmosphere is not well constrained by 
collider measurements (101). The various calculations differ by up to 2 orders of 
magnitude (18). Some of them are shown in Fig. 3 (59-61). The experiments 
LHCf and TOTEM, dedicated to measure particles emitted in the very forward 
region (102), will improve our understanding of hadronic models at laboratory 
energies of about 10 17 eV. 

A cubic-km NT measures of the order of kHz muon events, the so-called "sec- 
ond calibration beam." With such a high statistics, high precision verifications 
become possible; e.g., how well the detector details and the effects of the light 
propagation in the radiator are simulated. Nonetheless, hadronic interaction 
models used in the atmospheric shower simulations are affected by theoretical 
uncertainties, not only experimental ones. Some of these models are similar to 
those used to simulate atmospheric neutrinos. This represents a drawback, be- 
cause the two calibration beams should be used to understand the 'same' detec- 
tor independently of theoretical uncertainties. Another complication that affects 
shower models is the difficulty of investigating separately errors coming from im- 
perfect knowledge of high energy primary cosmic ray composition and hadronic 
model, and from the fact that simulations do not account for the seasonal vari- 
ation of the atmosphere. These simulations are very time consuming since rates 
of events are very high and high energy showers develop many particles in the 
atmosphere. Hadronic interaction models that foresee a larger K + production 
(such as SIBYLL 2.1) show better agreement than others (e.g. QGSJET II) in 
the high x = Secondary /-^primary region (103). The preference for such hadronic 
models is also indicated by the muon ratio measurement of MINOS (104) that 
shows a rise between 0.3 — 1 TeV consistent with an increasing contribution to 
the muon charge ratio from kaons. 

A neutrino telescope's high statistics measurement of the atmospheric neutrino 
flux above about 100 GeV is an important benchmark for all these calculations. A 
fit of the atmospheric muon neutrino flux [see Eq. (24)] was recently performed 
with AMANDA-II data to determine variations with respect to existing mod- 
els (105). The allowed region in the parameter space normalization vs spectral 
index was derived relative to the Bartol flux (98), using the angular distribution 
and an energy proxy for 5511 neutrino events collected in 1387 d. The best fit 
point indicates that, at 640 GeV, the data prefer a higher normalization by about 
10% and a slightly harder spectrum by about 0.056. Using a sample of 4492 upgo- 
ing neutrino induced muon events collected in 275.5 days, with a contamination 
of 5% from misreconstructed atmospheric muons, the IceCube Collaboration un- 
folded the atmospheric neutrino spectrum for the 22 string configuration (57). 
This spans a wider energy range up to PeV energies and is compared to some 
conventional (from n and K decays) and prompt neutrino models in Fig. 3. 
Systematic errors are so large at E u > 10 TeV that it is not yet possible to disen- 
tangle an eventual contribution from prompt neutrinos (originating in the decay 
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of charmed mesons and baryons) from an astrophysical component (60). 

Neutrino telescopes can collect a huge statistics of muons, of the order of bil- 
lions. Thus, they are sensitive to anisotropies of amplitudes < 1CT 3 . For ex- 
ample, in Fig. 10 we show the map in equatorial coordinates of about 5 billion 
atmospheric muons detected with the 22 strings of IceCube. These muons are 
produced by cosmic rays with a median energy of 14 TeV and have a median 
angular resolution of 3°, and show an anisotropy in their arrival direction of am- 
plitude (6.4 ± 0.2) x 10~ 4 . The observed anisotropy persists with high statistical 
significance, to energy in excess of 100 TeV with a smaller but measurable ampli- 
tude. This observation is not compatible with the Compton-Getting effect (109), 
which produces an energy-independent dipole anisotropy in the arrival direction 
of cosmic rays due to the relative motion of the solar system with respect to the 
rest frame of galactic cosmic rays (GCRs). This explanation was already excluded 
by Tibet Array (107) that in the absence of an anisotropy at 300 TeV concluded 
that GCRs co-rotate with the Galaxy. Other experiments observed anisotropies 
in the Northern hemisphere at the scale of 10° — 30° (108). Particularly Mi- 
lagro applied a cut that discriminates hadrons with respect to gammas. The 
observed anisotropy with IceCube in the opposite hemisphere persists at energies 
> 100 TeV, and hence makes it difficult to explain it as an effect of the helio- 
sphere or the heliospheric tail. Interestingly, it has been speculated that the effect 
may be due to the propagation of cosmic rays produced by nearby (~ 100 pc) 
SNRs and pulsars (110). It is noteworthy the presence of Geminga between the 
anisotropy regions in the Northern hemisphere. In the Southern hemisphere, the 
Vela pulsar is located near the highest intensity regions observed by IceCube. Of 
course such regions are so broad that we can only speculate about responsible 
local sources in the absence of a more precise calculation, including the intensity 
of galactic magnetic fields. This in turn could be correlated to observations at 
lower energies of an excess in the positrons to electron fraction with respect to 
secondary production models (111) that have also been interpreted as propaga- 
tion of cosmic rays from local sources such as SNRs and pulsars (112). It has 
also been speculated that the sites of the acceleration of cosmic rays may not be 
isolated, but also involve superbubbles (overlapping SNRs) (113). 

4.3 Searches for diffuse cosmic neutrino fluxes 

Diffuse fluxes from extra-galactic sources are more promising compared to single 
source fluxes because of the larger event rates, with a signal that may extend up 
to the highest energies where some of the extra-galactic sources (e.g., as AGNs 
and GRBs) can contribute. Nonetheless, these searches have to rely more on 
simulations than point source searches and hence are subject to larger systematic 
errors. The basic concept is that because of the harder spectrum the astro- 
physical signal should show up at high energies above the atmospheric neutrino 
background. However, the theoretical uncertainties on high energy atmospheric 
neutrino fluxes (see Sec. 4.2) and the experimental errors in the high energy re- 
gion where the statistics is low, limit the sensitivity and reliability of this search. 
For large statistics, the uncertainty in the prompt neutrino flux can be mitigated 
by triggering on tau neutrinos, which have the double bang characteristic signa- 



In Search of Extraterrestrial High Energy Neutrinos 



23 



ture in the PeV energy range, see Fig. 4. Specifically, the detection of a double 
bang event would confirm neutrino oscillations and at the same time establish 
the discovery of high energy cosmic neutrinos; this is because tau neutrino pro- 
duction in the atmosphere is negligible. Figures 11 (experiments measuring 1 
flavor) and 12 (experiments measuring all neutrino flavors) present a collection 
of experimental limits compared to representative models for classes of sources, 
such as AGNs, GRBs, and cosmogenic neutrinos. These upper limits can in turn 
be used to constrain the origin and nature of ultra-high energy cosmic rays. It is 
this that we now turn to study. 



4.4 Cosmic ray proton fraction 

Consider the working hypothesis introduced in Sec. 2 for optically thin sources, 
i.e., the relative number and energy of the neutrinos and neutrons depend only 
on kinematics, which implies approximate equipartition of the decaying pion's 
energy between the neutrinos and the electron, and the relative radiation density 
in the source. On average, each interaction will produce r\ neutrinos per neutron 
with relative energy e per neutrino, i.e. n = (N u ) / (N n ) and e = {E v )/{E n ). The 
neutrino emissivity of flavor i is then given by 

AR AR 

-j^-C Vt (z,E Wl ) = 1 ^£ n (z,E n ). (25) 

Assuming flavor universality as well as e ~ E u /E n ~ /S.E u //S.E n and n ~ 
Ngjiv/Nn one arrives at the neutrino source luminosity (per co-moving volume) 

C^ v (z,E v ) ~ ^Cn(z,E v /e). (26) 

For the (hypothetical) source where pion production proceeds exclusively via res- 
onant p7 — > A + scattering, the neutrino multiplicity and the relative energy are 
fixed: r] = 3 and e = 0.07. Note, that the relation (26) derived for optically thin 
sources can be regarded as a lower limit on the neutrino luminosity as long as 
energy-loss processes in the source are negligible (122). Therefore, this conserva- 
tive expectation can be translated into an upper limit on the extra-galactic proton 
fraction in ultra-high energy cosmic rays, exploiting, e.g., the experimental up- 
per limit on the diffuse neutrino flux from IceCube or the preliminary AMANDA 
data (119). 

For this procedure one introduces test functions of the neutron source luminos- 
ity of the form C^ st (0,E) = Cq (E / E m&yi )~ l exp(— E/E max ), with an exponential 
energy cut-off E mSLX that varies between 10 8 GeV and 10 12 GeV with a logarith- 
mic step-size of Alog 10 E = 0.25 (122). The source function per unit redshift is 
given in Eq. (19). The source luminosity per co-moving volume is assumed to 
follow that of the star formation rate (SFR): Hsfr(^) = (1 + Z ) 3A 7 f° r z < 1, 
Wsfr(^) = Aq(l + z)-°- 3 , for 1 < z < 4, and H SFR (z) = AqTV^l + z)~ 3 - 5 , for 
z > 4, with appropriate normalization factors, N\ = 2 3 ' 7 and N4 = 5 3 ' 2 (123). 
Each neutron test luminosity is related to a neutrino luminosity by the ratio (26). 
After propagation, the accumulated contribution of extra-galactic and cosmogenic 
neutrinos is normalized to the neutrino limit. The results are shown in Fig. 13. 
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The AMANDA preliminary limit on the diffuse neutrino flux (119) severely con- 
strains the extragalactic proton contribution in cosmic rays at energies below 
10 9 GeV. Following the same line of thought, one can verify that if the evolution 
of cosmic ray sources follows that of active galactic nuclei (124), the upper limit 
on the diffuse flux of tau neutrinos from the Pierre Auger Observatory (117) 
marginally constrains the proton fraction at the end of the spectrum (125). 

4.5 Dark matter searches 

Over the past few years, a flood of high-quality data (126) from the Super- 
nova Cosmology Project, the Supernova Search Team, the Wilkinson Microwave 
Anisotropy Probe (WMAP), the Two Degree Field Galaxy Redshift Survey (2dF- 
GRS), and the Sloan Digital Sky Survey (SDSS) pin down cosmological parame- 
ters to percent level precision, establishing a new paradigm of cosmology A sur- 
prisingly good fit to the data is provided by a simple geometrically flat Friedman- 
Robertson- Walker universe, in which 30% of the energy density is in the form of 
non-relativistic matter (f2 m = 0.30±0.04) and 70% in the form of a new, unknown 
dark energy component with strongly negative pressure (Oa = 0.70 ± 0.04). The 
matter budget has only 3 free parameters: the Hubble parameter h = 0.70+QQ3, 
the matter density fi m /i 2 = 0.138 ± 0.012, and the density in baryons, tl^h 2 = 
0.0230^q'qq];2- 5 This implies that the structure of the universe is dictated by the 
physics of as-yet-undiscovered cold dark matter (Ocdm^ 2 = 0.115 ± 0.012) and 
the galaxies we see today are the remnants of relatively small overdensities in the 
nearly uniform distribution of matter in the very early universe. 

The simplest model for cold dark matter consists of WIMPs - weakly inter- 
acting massive particles (128). The most popular CDM candidates are (i) the 
lightest supersymmetric particle (such as the well motivated neutralino %) that 
is stabilized in i?-parity conserving models of supersymmetry (SUSY) (129) and 
(ii) the lightest Kaluza-Klein state appearing in models of universal extra dimen- 
sions (130). Generic WIMPs were once in thermal equilibrium, but decoupled 
while strongly non-relativistic. The relic abundance of SUSY WIMPs can be 
found by integrating the Boltzman equation (131), 

d ^ +2> Hn = -{av){n 2 -nl q ), (27) 

where n is the present number density of SUSY WIMPs, H is the expansion rate 
of the universe at temperature T, n eq = g (m x T/27r) 3//2 e~ m ^l T is the equilibrium 
number density, g is the number of internal degrees of freedom of the WIMP, (<jv) 
is the thermally averaged annihilation cross section, and m x is the neutralino 
mass. Note that in the very early universe, when n ~ n cq , the right hand side 
of Eq. (27) is small and the evolution of the density is dominated by Hubble 
expansion. As the temperature falls below m x , however, the equilibrium number 
density becomes suppressed and the annihilation rate increases, rapidly reducing 
the number density. Finally, when the number density falls enough, the rate 
of depletion due to expansion becomes greater than the annihilation rate, H > 

5 The latter is consistent with the estimate from Big Bang nucleosynthesis, based on mea- 
surements of deuterium in high redshift absorption systems, flbh 2 — 0.020 ± 0.002 (127). 
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n{crv), and the neutralinos freeze out of thermal equilibrium. The freeze-out 
temperature Tf depends logarithmically upon (crv), but for models with TeV 
scale SUSY breaking, one finds that T { /m x ~ 0.05. SUSY WIMPs are, by far, 
the favored candidate for CDM, because for masses of order 100 GeV to 10 TeV 
a present density of Q x h 2 ~ 0.1, comes out fairly natural (132). 

Many approaches have been developed to attempt to detect dark matter. Such 
endeavors include direct detection experiments which hope to observe the scat- 
tering of dark matter particles with the target material of the detector (133) and 
indirect detection experiments which are designed to search for the products of 
WIMP annihilation into gamma-rays, anti-matter, and neutrinos (134). 

The annihilation of WIMPs into photons typically proceeds through a com- 
plex set of processes, and the final spectrum actually contains a detailed imprint 
of WIMP annihilation that can, in principle, reveal features such as the WIMP 
spin and/or other particles in the dark sector (135). On the one hand, tree-level 
annihilation of WIMPs into quarks or leptons (or heavier states which decay 
into them) yields a continuum of photon energies, with an upper cutoff at ap- 
proximately the WIMP mass. On the other hand, loop-level annihilation into a 
photon and X results in a line at energy E^ = m x [l — M|/(4mp, where Mj is 
the X-boson in the final state (in SUSY models, X is either another photon or a 
Z). The line emission typically has smaller magnitude than continuum emission, 
but it provides a clean signal that helps discriminate against backgrounds. The 
differential flux of photons arising from CDM annihilation observed in a given 
direction making an angle ip with the direction of the galactic center (GC) is 
given by (135) 

W.W4^|£M^<«, (28, 

where J = (1 /AS7) J A ^ J{^) dQ, denotes the average of J over the solid angle Af2 
(corresponding to the angular resolution of the instrument) normalized to the 
local density: J(ip) = (DqPq) -1 f^ Q p 2 [r(£,ip)]d£; the coordinate £ runs along 
the line of sight, which in turn makes an angle ip with respect to the direction of 
the GC ( i.e., r 2 = I 2 + Dq — 2£D & cos ip)', the subindex / denotes the annihi- 
lation channels with one or more photons in the final state and dNf/dE^ is the 
(normalized) photon spectrum per annihilation; and p(x), = 0.3 GeV/cm 3 , 
and Dq ~ 8.5 kpc respectively denote the CDM density at a generic location x 
with respect to the GC, its value at the solar system location, and the distance 
of the Sun from the GC. Indeed, the GC has long been considered to be among 
the most promising targets for detection of CDM annihilation, particularly if the 
halo profile of the Milky Way is cusped in its inner volume (136). This has been 
complicated, however, because H.E.S.S. observations disfavored CDM emission 
from the GC, and revealed instead a VHE source (HESS J1745-290) (137). To 
discover CDM photon emission the contribution of the many point sources near 
GC (138) needs to be subtracted. IACTs or satellite and ballon-borne experi- 
ments also look for gammas from high CDM density regions relatively close to 
the Earth (> 100 kpc) such as dwarf spheroidal satellite galaxies of the Milky 
Way, clusters of galaxies, and intermediate mass black holes (85). 
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Neutrino telescopes indirectly search for the presence of dark matter by taking 
advantage of the Sun's ability to capture large numbers of WIMPs over time. 
Over billions of years, a sufficiently large number of WIMPs can accumulate in the 
Sun's core to allow for their efficient annihilation. Such annihilations produce a 
wide range of particles, most of which are quickly absorbed into the solar medium. 
Neutrinos, on the other hand, may escape the Sun and be detected in experiments 
on the Earth. Specifically, the time evolution of the WIMP population in the Sun 
is controlled by 

N = C e - A e N 2 , (29) 

where 

C -SxlO^s- 1 ( Aocal ^ f 270 km/* W geff \ ATeV\ 2 
Ce ~ 3Xl ° S U-3 GeV/cm3j V^^Tj [lO^ph) ) (30) 

is the Sun's rate of capture (139), pi oca i is the local dark matter density, v\ oca x 
is the local rms velocity of halo dark matter particles, Aq = (crv)/V e e, V e e = 
1.8 x 10 26 (TeV/m x ) 3 / 2 cm 3 is the effective volume of the core Sun determined by 
matching the core temperature with the gravitational potential energy of a single 
WIMP at the core radius (140), and A G is related to the annihilation rate by 

V= l -A Q N 2 . (31) 

The rate at which WIMPs are captured in the Sun depends on the nature of 
the interaction the WIMP undergoes with nucleons in the Sun. These elastic 
scattering processes are often broken into two classifications: spin dependent 
interactions in which cross sections, ctsd> increase with the spin of the target 
nuclei, and spin independent interactions in which cross sections, ctsi, increase 
as the square of the total number of nucleons in the target. Thus, the effective 
WIMP-on-nucleus (hydrogen and helium) elastic scattering cross section has two 
contributions: a e s = o"sd + cid- The present WIMP annihilation rate is found by 
solving Eq. (29) for the annihilation rate at a given time 

r = ^ tanh 2 (-£) , (32) 

where r eq = \j \JCq Aq is the time scale required to reach equilibrium between 
capture and annihilation. When r eq becomes comparable or larger than the age 
of the solar system (i Q ~ 4.5 billion years), the system has not yet reached equi- 
librium and the annihilation rate is strongly suppressed r \ Cq(^Cq A & t G ) 2 , 
whereas when (Cq Aq) i Q 3> 1 the annihilation rate is saturated at T = \Cq. 

As they annihilate, WIMPs can generate neutrinos through a wide range of 
channels. Annihilations to heavy quarks, tau leptons, gauge bosons and Higgs 
bosons can all generate neutrinos in the subsequent decay (22). The total flux of 
neutrinos emitted by the Sun due to WIMP annihilition is then 

= ^ q £(^). e -^ /150GeV , (33) 



dF^ 



dE v 
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where F^ q is the non-equilibrium suppression factor (~ 1 for capture-annihilation 
equilibrium) and (dN Ufi /dE u )j is the flux produced by the j decay channel per 
WIMP annihilation. Note that neutrinos produced near the center of the Sun 
interact with the solar medium, yielding a depletion of the emission spectrum by 
a factor ~ e -E v /i50 GcV (22) . Finally, the flux reaching the Earth is, 



1 dF M 

k"» = ® 

Aitd 2 dE u 



(34) 



where d « 1.5 x 10 8 km is the Earth-Sun distance. All in all, the prospects 
for neutrino-detection-experiments detecting dark matter critically depend on 
the capture rate of WIMPs in the Sun, which in turn depends on the elastic 
scattering cross section of these particles. 

While gamma-ray experiments set limits on (av) or on the gamma-ray flux, 
neutrino telescopes measure the muon induced flux by neutrinos and hence con- 
strain this flux as a function of m x . A collection of various limits on the muon flux 
is shown in Fig. 14. The conversion from muon flux to cross section limits is not 
model independent; in Fig. 15 it was done using DarkSUSY (150). Equilibrium 
is assumed between capture and annihilation rates in the Sun, so that the annihi- 
lation rate is proportional to the spin-dependent and independent cross sections. 
A limit on o"sd is found by setting to zero the spin-independent cross section. 
This procedure is indispensable to show complementarity between searches and 
to combine results from different techniques. As an illustration, in Fig. 15 we 
compare the limits on ctsd with those obtained from direct detection experiments 
including CDMS (141), COUPP and KIMS (149). 

The described indirect searches for CDM are complementary approaches to 
the direct searches and the possible direct production at LHC. Though indirect 
detection channels are subject to large astrophysical uncertainties, they explore 
different regions of the parameter space and a very wide range of candidate par- 
ticle masses. 



5 CONCLUSIONS AND FUTURE PROSPECTS 

The hunt for neutrinos has just begun with suitable size detectors for the expected 
small event rates. Given their dimensions these experiments are very challenging 
to build. IceCube, the first cubic-kilometer detector, is expected to be completed 
on schedule in 2011. If the predictions from the measured gamma-ray fluxes are 
indeed correct, then the first detection of sources, such as the "Milagro Pevatrons" 
with hard spectra up to a few tens of TeV, may become possible after a few years 
of data taking. Neutrino telescopes may broaden our reach to the unknown like 
many other astronomical instruments which have discovered unexpected signals 
from the Universe; NTs can access sources that are opaque in photons. The 
Mediterranean community is operating ANTARES, the first complete underwa- 
ter detector, at the scale of the IceCube's precursor AMANDA-II. Other R&D 
programs are ongoing to build a cubic-km detector in the Mediterranean sen- 
sitive to the GC. In Italy, a ~ 100 km electro-optical cable has been installed 
off-shore Capo Passero (at about 3.5 km below sea surface) by the NEMO Col- 
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laboration, in July 2007. A prototype tower of the proposed cubic-km array will 
be connected in the near future. ANTARES, NEMO and NESTOR joined their 
efforts forming the Consortium KM3NeT to design the cubic-km detector in the 
Mediterranean. High quantum efficiency PMTs are being studied for this cubic- 
km detector in the Mediterranean. Most of the proposals aim at increasing the 
photocathode area per OM and the isotropy of light detection, minimizing the 
transit time spread, maximizing the collection efficiency and the peak to valley 
ratio for a single photoelectron signal (connected to the capability of counting 
photoelectrons in one event), and minimizing costs and power consumption. 

Building a detector a factor of 10 more sensitive than IceCube does not seem to 
be an easy task to achieve, neither in the ice nor in the sea. A new detector in the 
sea, with similar cost to IceCube, could not achieve such a goal with the current 
technology, despite the fact that in the sea the angular resolution for point source 
searches can be better by a factor of 2-3 compared to IceCube. It is also hard to 
obtain a uniform increment in area at all energies, while avoiding cost divergence. 
A significant increase in the area at very large energies implies larger separations 
between instrumented lines. Hence, this penalizes the 1—100 TeV region, that is 
an important region for galactic sources. 

Projects are underway to build detectors of the order of 100 times IceCube 
using techniques that exploit the advantage of attenuation lengths about a factor 
of 10 greater than the optical one in the same medium. Particularly interest- 
ing is the radio technique in ice. The acoustic technique in sea water or ice 
seems also encouraging. None of these techniques, however, is background free. 
Moreover, they have the disadvantage of higher energy thresholds compared to 
the NTs discussed in this review. Clearly the higher thresholds reduces the de- 
gree of overlap with the optical technique, preventing effective cross-calibrations 
that would help the understanding of the background. Moreover, thresholds of 
the order of 10 18 eV are suitable for extragalactic neutrinos but not for galactic 
sources. Detectors combining the radio and acoustic techniques are also desir- 
able. However, for cross-calibrations, the radio technique seems currently more 
promising due to the low cost of antennas, the lower power consumption than 
acoustic devices, and the larger attenuation length in the ice. 

In summary, after 20 years of careful work by many research groups around 
the world we are in possession of a tantalizing body of data, more than sufficient 
to stimulate our curiosity but not yet sufficient to pin down high energy neutrino 
sources. The upcoming high quality observations from the "giant-aperture" neu- 
trino telescopes under construction will generate a data sample of unprecedented 
size, ushering us to a golden age of neutrino astronomy. 
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Figure 1: The solid lines indicate the range for the cosmogenic neutrino spectrum 
for a = 1.4 — 3.0 and E max /Z = 10 21 — 10 22 eV, for pure iron (up) and iron and 
proton (down). In each frame, we show for comparison as a dashed curve the 
prediction for an all-proton spectrum with a = 2.2 and £ max = 10 22 eV (16). 



40 Anchordoqui Sz Montaruli 



is 10 



ANTARES 1 line 
AMANDA-II 2003 
DUMAND 1990 
Baikal NT-36 1993 
NESTOR prototype 89-92 

NESTOR prototype 2005 
Sinegovskaya et al 2001 E >2Q GeV 
Sinegovskaya et al 2001 e" >1 GeV 




4 5 6 7 

Depth [km.w.e] 







Figure 2: Atmospheric muon vertical intensity measured by under- water and ice 
arrays as a function of depth compared to calculations (52). No correction is 
applied for the slightly different muon propagation properties between water and 
ice. Data are from: ANTARES and AMANDA-II (53), Baikal, NT-36 (54), the 
DUMAND prototype string (55), and NESTOR (for 2 different prototypes) (56). 
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Figure 3: Preliminary atmospheric unfolded neutrino spectrum measured with 
22 strings of IceCube (57) compared with predictions for z/„ + (The factor 

1/3 accounts for the lower CC interaction cross section of anti-neutrinos up to 
about 10 5 GeV compared to neutrinos.) Data are accompanied by large error 
bars. A preliminary evaluation of systematic errors, which are mainly due to the 
ice purity level dependence with depth, is included in the error bars. The energy 
bin size was selected to match the energy resolution of the experiment. (Since the 
energy resolution is 0.3 on the log 10 (-E / i), there are 3 bins per decade.) The size 
of the last 2 bins is larger, and determined by the desire to accommodate bins 
with vanishing contents. Conventional neutrino fluxes are calculated as in (58), 
prompt models are calculated in the framework of perturbative-QCD in (60) 
(std stands for optimal parameters, whereas min and max indicate the range of 
variation of the model parameters) and in (59) for different structure functions; 
the Quark-Gluon String model (QGSM) (61) is also shown. 




Figure 4: Simulated events in the IceCube detector, visualized using the IceCube 
event display, showing the 3 typical topologies discussed in Sec 3. The shading 
represents the time sequence of the hits. The size of the dots corresponds to the 
number of photoelectrons detected by the individual photomultipliers. From left 
to right: a muon event of 100 TeV, a cascade event induced by a 100 TeV u e , 
and a double bang event induced by a 16 PeV u T . 
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Figure 5: The IceCube Observtory, including the deep ice array, Ice Top, 
AMANDA, and Deep Core. For comparing sizes the image of the Eiffel tower is 
also shown. 
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Figure 6: A pictorial view of the ANTARES detector including an insert of a 
storey. 




Figure 7: Effective area for muon neutrinos, integrated over the lower hemisphere, 
as a function of the neutrino energy. The selected cuts ensure a pointing capability 
suitable for point source searches. The area for 40 strings of IceCube is obtained 
with the requirement that tracks are reconstructed inside 2° from a source. The 
full IceCube area is obtained using the same cuts, and is therefore preliminary. 
The AMANDA-II area is from the final analysis in Ref. (80). The ANTARES 
effective area is obtained with preliminary selection cuts based on simulation (81). 
The KM3NeT area is for a possible configuration still under debate (82). 
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Figure 8: Cumulative distribution of reconstructed muon events inside a cone of 
semi-angle (Point Spread Function) from the neutrino direction. The dif- 
ferent curves indicate growing configurations of IceCube for two neutrino energy 
bins and the point-source analysis cuts (83). The IceCube curves are preliminary, 
because they are based on cuts optimized for the 40 string analysis. 




Figure 9: Muon event skymap (z-axis color code is pre-trial significance) for six 
months of data taken with the 40-string IceCube array, from July 2008 through 
December 2008. Of the 17777 black dots on the skymap, 6797 are upgoing neu- 
trino candidates (zenith > 90°) from the northern hemisphere. The color shading 
indicates the significance of the data and the curved black line is the galactic 
plane (83). 
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Figure 10: Upper panel: Celestial cosmic ray intensity map for Tibet AS7 data 
taken from 2001-5. The vertical color bin width for the relative intensity, that 
is calculated within each declination band, is 2.5 x 10~ 4 (107). The highest 
and lowest intensities are observed in the region of the tail-in and loss-cone in the 
direction of the heliotail, and on the right the excess is in the direction of Cygnus. 
Lower panel: Downward going muon skymap (expressed as relative intensity 
mini tier of events in an angular bin in right ascension and declination divided by 
the total number of events in the declination band) observed by 22 strings of the 
IceCube neutrino telescope with data taken between April 2006-7 (106). 
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Figure 11: The blue triangles indicate the atmospheric {y^ + £v) flux averaged 
over the lower hemisphere as measured by 22 strings of IceCube (57)), and the 
falling thick solid line indicates AMANDA-II measurements (105)). For com- 
parison, predictions of the atmospheric neutrino flux (conventional + prompt) 
are indicated by the falling thin solid lines (58,60,61,98). The horizontal lines 
represent different 90% CL limits and sensitivities on E~ 2 fluxes; from top to 
bottom: ANTARES estimated sensitivity for 1 yr (114), AMANDA-II muon 
neutrino limit for 804 days (115), IceCube estimated sensitivity for 1 yr (116). 
For comparison, also shown are (%) the Waxman and Bahcall upper limit cor- 
rected for oscillations (42), (ii) neutrino flux predictions from AGNs (12) and 
prompt and precursor emission from GRBs (13), and (Hi) various extimates of 
the nearly guaranteed cosmogenic flux (14, 15). The ultra-high energy upper limit 
on the tau neutrino flux (differential and oc E~ 2 ) reported by the Pierre Auger 
Collaboration is also shown (117). 



In Search of Extraterrestrial High Energy Neutrinos 47 



1.1 (T 



Atmospheric v„+v e max 
Atmospheric v tl +v c min 
AMANDA-I1 1001 d prelim 
IceCube 22 strings EHE 242 d 
AMANDA-II UHE 507 d prelim - 
Allard et al, 2006, protons 
Allard et al, 2006, protons 
ESS v„ + v e , 2001 
HBL Mucke et all, 2003 x 3 



Rice 854 d 
HiRes v T +v e 
Baikal 1038 d 
ANITA-2008 35 d 
Waxman and Bahcall 1998 x 3/2 
Razzaque et al GRB Progenitor x3/2 
Razzaque et al Prompt GRB x3/2 
Blazars Stecker 2005 x 3 




10 12 14 

log io E v [GeV] 



Figure 12: Same as in Fig. 11 but for experiments measuring more than one 
flavor, e.g. measuring cascades induced by all flavor neutrinos. The horizontal 
solid lines are 90%CL limits for a flux cx E~ 2 ; from top to bottom: AMANDA-II 
all-flavor cascade limit for 1001 days (118), Baikal cascade limit for 1038 days (69), 
AMANDA-II muon neutrino limit for 804 days (115), all flavor ultra-high energy 
limit (preliminary) for 507 days of AMANDA-II (119). The differential in energy 
upper limits (90% CL) are for the 22 strings of IceCube (120) and several ultra- 
high energy neutrino detection experiments, including those with radio detectors 
that look for the coherent Cherenkov radiation produced by the excess of electrons 
in neutrino induced showers in dense media (e.g., ice and the lunar regolith). For 
further details see (85). 
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Figure 13: Integrated and differential upper limits on the proton contribution in 
ultrahigh energy cosmic rays derived from AMANDA bound on diffuse neutri- 
nos (119). (After Ref. (122).) 
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Figure 14: Upper limits (90%CL) on the muon fluxes from neutralino annihila- 
tion in the Sun as a function of neutralino mass for soft (bb) and hard (W + Ty~) 
channels (85, 147). The lighter hatched region indicates the SUSY parameter 
space compatible with direct detection limits on the spin independent cross sec- 
tion from CDMS (141) and XENON10 (142). The darker hatched area rep- 
resents the projected sensitivity for the same region of the parameter space, 
assuming that direct detection limits are a 100 times better. Experimental lim- 
its, that include a correction for detector threshold within the common assump- 
tion of E V) fhi = 1 GeV, are shown for Super-Kamiokande (143), MACRO (144), 
AMANDA-II (145), Baikal (146) and IceCube 22 strings for the hard and soft 
channels (147). For comparison, the projected sensitivity of ANTARES after 5 yr 
of data taken is also shown (148). 
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Figure 15: 90% CL upper limits from AMANDA-II and 22 strings of IceCube 
on the spin dependent cross section for soft and hard channels (147). Hatched 
regions follow the same conventions of Fig. 14. For comparison, we also show 
direct search limits on cj S d from CDMS (141), COUPP (149) and KIMS. 



